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Abstract 

The content of this Quarterly Report can be summarized as follows: 

a) A brief study has been carried out on a retrieval rate control law with 
no angular feedback in order to investigate the dynamic response of the system. 

b) The initial conditions for the computer code which simulates the 
satellite's rotational dynamics have been extended to the case of a generic 
orbit . 

c) The model of the satellite thrusters has been modified in order to 
simulate a pulsed thrust (actual case) , by making the SKYHOOK integrator 
suitable for dealing with delta functions without loosing computational 
efficiency. 

d) Tether breaks have been simulated with the high resolution computer code 
SLACK3. Shuttle's maneuvers in order to avoid the recoiling tether have been 
tested . 

e) The electric potential (in vacuo) around a severed conductive tether 
with insulator, in the case of a tether breakage at 20 km from the Shuttle, has 
been computed. 

f) Preliminary evaluation of the electrodynamic hazards due to the breakage 
of the TSS elactrodynamlc tether in a plasma has been carried out. 
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1.0 INTRODUCTION 

This is the fourth Quarterly Report submitted by SAO under contract 
NAS8-36160, "The Investigation of Tethered Satellite System Dynamics," Dr. 
Enrico Lorenzinl, PI, and covers the period from 15 May 1985 through 14 August 
1985. 


2.0 TECHNICAL ACTIVITY DURING REPORTING PERIOD AND PROGRAM STATUS 

2.1 Retrieval Using Rate Control 

There are two basic types of retrieval control laws - tension control and 
rate control. The results reported in Quarterly Report 03 used tension control. 
Tension control laws are easily integrated into programs such as DUMBEL and 
SKYHOOK which integrate the motion in inertial coordinates given the forces on 
each mass point. In computer programs where the length of the wire is one of 
the integration variables, rate control laws can be easily Implemented by using 
a first order differential equation for the length with the rate of change of 
length given by the control law. Such a system has a potential disadvantage in 
that there is no guarantee of positive tension in the wire since the radial 
acceleration and the tension are not involved in the integration. Rate control 
laws can be implemented in programs such as DUMBEL and SKYHOOK by including the 
natural length of the wire among the variables to be integrated, The rate of 
change of the length is specified by the control law. The tension in the wire 
is computed from the difference between the actual length and the natural 
length. In this system one can set the tension to zero if the actual length is 
less than the natural length so that there is no possibility of having negative 


tension in the simulation. This method of simulation has a high fidelity to the 
actual physical situation where the reel motor pulls in wire, increasing the 
tension and accelerating the end mass . 

The DUMBEL program already has a facility for integrating the natural 
length of the wire. Under a previous contract a model of a damper was added to 
the program in which the length of wire released by the damper reel is included 
as an integration variable. The program has been modified to use this variable 
for rate control in the retrieval mode of the program. For a perfect retrieval 
at a constant in-plane angle 0 C , the rate of change of wire length is 

t - - at (1) 

where 

3 

a — — w 0 cos 0 C sin 0 C (2) 

z 

and w 0 = the orbital angular velocity. 

A rate control retrieval mode has been added to the DUMBEL program using 
equation (1) to specify the rate of change of the natural length of the wire. A 
total of thirteen variables are integrated in the version without rotational 
dynamics. Six variables are required for the position and velocity of each mass 
and one for the variable giving the change in the natural length. 

A test run has been done using equation (1) for retrieval in the DUMBEL 
program without rotation. The integration was started with equilibrium initial 
conditions consisting of a 7 m/sec retrieval velocity and a 12.1846 deg in-plane 
angle with 0 o — 12.1846 deg. The subsatellite was deployed downward on a 20 km 
tether. Figure 2.1.1 shows a plot of the in-plane angle vs. time for the first 
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15,000 seconds. The initial angle of 167.8152 is 12.1848 deg from the vortical 
position at 180 degrees. At 13,600 sec the subsatellite was 20.2179 degrees 
from the vertical. At 15,000 sec it is about to pass 'through the vertical at 
180 degrees, As can be soon in the figure the ln-plano libratlons are 
increasing in amplitude during the retrieval. The wire length at 15,000 seconds 
was 104 meters. The run was continued to 20,000 seconds and a wire length of 18 
meters. The maximum in-plane deviation from the nominal retrieval angle was 
about 50 degrees. Although the run did not go unstable in the sense of having 
the subsatellite wrap-around the Shuttle, it is clear that this control law does 
not provide any damping to prevent the buildup of in-plane libratlons. 
Non-equilibrium Initial conditions would presumably show more rapid development 
of libratlons, A damping term could be added to the rate control law using 
angle feedback. However the necessary information is not expected to be 
available to the reel control system on the first mission. 


2.2 Initial Conditions For Inclined Orbits 

The DUMBEL program has a facility for generating equilibrium initial 
conditions with the tether displaced from the vertical in the in-plane or 
out-of-plane directions, The original software was written in a simplified 
manner which could only be used for equatorial orbits. This facility has been 
updated by including an additional rotation to account for the orbital 
inclination. This facility will be needed to do simulations to study the 
effects of earth oblateness, etc, A short test run was done for a 28° 

inclination orbit with a 10 degree initial in-plane displacement and no 
out-of-plane displacement. Plots of the llbratlon angles vs, time showed the 
proper behavior. The software changes were made in subroutine LAUNCH which 


reads the libration angles at?d modifies the state vector, 


The format of the output on DUMBEL has been changed to label each variable 
at each output point. The old format with a single heading at the beginning of 
the run was confusing to road on later pages. 


2.3 Use Of Thrusters For Controlling Librations 
2.3,1 Software Implementation - 

The TSS is expected to have thrusters in the in-plane and out-of-plane 
directions which can be fired manually by commands from the Shuttle. In 
addition there will be an in-line thruster for use in maintaining wire tension 
during the final stages of retrieval. The thrust level for the transverse 
thrusters is 2.5 newtons and for the in-line thruster 2.0 newtons. The operator 
can enter a command for the thruster to fire at a certain time for a given 
duration not less than about 50 milliseconds. In contrast to the tension 
control algorithm which can operate in an automatic mode using data on length 
and length rate available to the deployer, the control of librations using 
thrusters must be done manually since angle Information available in the cockpit 
is not transmitted to the deployer. Studies have been done to develop a 
procedure for using the thrusters to eliminate librations during retrieval. 

The work done by a thruster is proportional to the velocity of the 
subsatellite at the time of firing. Therefore the most efficient technique is 
to fire the thruster near the time when the system is swinging through the 
vertical position and is at maximum swing velocity, The periods of the in-plane 
and out-of-plane librations are about 3206 and 2776 seconds respectively at 400 
km altitude. The times for one quarter of a cycle are 801 and 694 seconds 


Pago 9 


respectively. In order to keep the firing near the maximum velocity and upper 
limit of 100 seconds has boon used for the firing time. 

If the velocity is low enough, the swinging motion can bo stopped with one 
firing. The firing interval required can bo computed, at least approximately, 
from the amplitude and period of the oscillation. Suppose the oscillation is 
given by the equation s 


x «= A sin (wt) ( 3 ) 

where x is the displacement of the subsatellite from the local vertical, A is 
the maximum displacement and w is the frequency of the libratlon. For the 
in-plane component, the frequency is y'T" and for tho out-of -plane 2 w 0 where 
w 0 is the orbital frequency, The velocity is 

x >= A to cos (wt) 


The maximum velocity is 


v = A (J 


The time required to stop the motion is 


(4) 


T = v/a (5) 

where a is the deceleration. If the mass of the subsatellite is m and the 
thruster force is F then the acceleration is 
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a a F/ra (6) 
Combining aquations 4, 5, and 6 gives 

T a m A w/F (7) 

The computar implementation of the firing commands must take into account 
the way the numerical integrator works. The integrator uses variable step size 
to maintain the integration accuracy requested by the user. The user must 
supply a subroutine that computes the forces on each mass as a function of the 
oxtrapolatod state vector computed by the integrator. Before taking a step, the 
integrator will try various step sizes to see how large a step can be taken and 
still maintain the required accuracy. The times given to the subroutine 
supplied by the user are therefore not monotonic and the subroutine must be 
prepared to deal with random input times. The stopsizo varies greatly depending 
on the nature of the function being integrated. At a discontinuity the stopsize 
may be as small at 10‘ 7 seconds. During periods when the behavior is smooth the 
stepeize may bo 25 seconds or more. 

A subroutine has been written to generate firing intervals for the 
thrusters. The subroutine operates in three modes. In the first mode the 
subroutine tests the libration amplitude to see if it has exceeded the deadband 
specified by a parameter. If the amplitude of the libration has exceeded the 
deadband a flag is set to search for the next maximum in the libration angle. 
At each call, the amplitude is compared to that on the previous call and the 
time and amplitude saved if the value is larger. If it has been more than 1/8 
of the libration period since the amplitude has increased the subroutine sets a 
flag indicating that a maximum has been reached and the program goes to the next 
mode, In the third mode the program searches for a zero crossing. This is done 
by lookin'? for a change in sign of the amplitude of the libration. Because of 


Pago 11 


the way the integrator works tho time of the zero crossing may La as much as 25 
seconds away from tho first time a change in sign is detected. In this 
preliminary version of the software the difference has been ignored. Tho 

subroutine simply takes the first time there is a sign change and sets this os 
the on time for the thruster, Tho firing interval is computed from equation 7 
using the maximum amplitude A determined in mode 2 of the subroutine. If the 
time T is greater than 100 seconds tho firing interval T is sot to 100 soconds. 
This system does not center the firing time around the zero crossing and the on 
time may be off by an amount not exceeding the integrator stepsize. 

2.3,2 Retrieval Simulation Using Thrusters - 

As a first tost of tho thruster control software the tether was given an 
initial in-plane displacement of 30 degrees. The ln-plane thrusters wore flrod 
each time the subsatelllte crossed the local vortical. The simulation was run 
for 10,400 soconds with no retrieval. The subsatellite was deployed downward on 
a 20 km tether. The thruster fired for 100 second intervals beginning at 859, 
2568, 4270, and 5968 seconds. The firings were effective in reducing tho 

amplitude of the libratlons. The approximate times and amplitudes of tho 
maximum libration angles are listed in Table 2.3.1. 

Table 2.3.1 

Maximum libration angles vs. time. 


T (sec) 

Amplitude 

0 . 

+ 

10.000 

1700. 

- 

9.768 

3400. 

+ 

9.535 

5100. 

- 

9.305 

6800. 

+ 

9.077 

8500. 

- 

8.850 


In the second tost run, the tether was given an initial out-of-plane 
displacement of 10°, The out-of-plane thrusters were fired for 100 second 
intervals each time the out*of-plane displacement went through zero. The 
amplitude of the out-of-plano displacement decreased at each thruster firing. 
However, the decrease was not sufficient to keep the out-of-plano angle from 
increasing during the retrieval, At around 8900 seconds, the radial component 
went from negative (downward deployment) to positive (upward deployment) . The 
tether length at that time was about one kilometer. From that point on, the 
subsatel lite was being retrieved from above. The run Was terminated at 19,500 
seconds at a tether length of about 47 meters. At the beginning of the run, the 
firing time computed from equation 7 was 1729 seconds. The computed firing time 
decreased after each thruster firing but by more than 100 seconds the largest 
decrease being about 337 seconds at around the time the system flipped to upward 
deployment. Since the amplitude of the llbratlon became very largo, the tension 
control law probably provided significant damping of the out-of-plane component 
since the coupling between the radial and out-of-plane components increases with 
amplitude. Toward the end of the run the firing interval dropped below the 
integration stepsize of 25 seconds, and the program started skipping over the 
firings, so that there was no further decrease in the out-of-plane component. 

The plot of the in-plane angle in the previous runs was not satisfactory in 
the case of downward deployment. The angle is 180° near the vertical, but the 
usual normalization is from -180° to +180° so that there is a change of 360° as 
the in-plano angle crosses the local vertical. The normalization was changed to 
be from -90° to +90° when the cosine of the angle is positive and 90° to 270° 
when the cosine of the angle is negative. 

In the next test run, the system was given ann initial out-of-plane 
displacement of 3° and the out-of-plane thrusters were used during retrieval. 
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The thruster wore fired for 100 seconds at each zero crossing until tho firing 
time required to stop the litaration went below 100 seconds. The firing times 
and the firing interval required to stop the libration are listed in Table 2.3.2 
for each zero crossing. 


Table 2.3,2 

Start Time for Each Thruster Firing and 
Firing Interval Required to Stop *,he Libration 


Thrustor Start Time 

660,1 

1997.9 

3432.5 

4821.3 

6208.3 
7479.1 


Computed Firing Duration 

521.1 

425.8 

345.3 

234.4 

136.5 
35.3 


The computed firing duration decreases by roughly 100 seconds on each zero 
crossing. The run was terminated at 7513.7 seconds because the integration 
stepslze became prohibitively small, In an attempt to understand the source of 
the problem, the run was repeated with output Intervals of 1 second starting at 
7480 seconds and continuing to the time at which the integration failed, The 
output interval on the original run was 100 seconds which did not give enough 
time resolution to study what was happening. Plotting the out-of -plane 
displacement vs. time showed that the libration velocity was going to zero at 
the point where the problem occurred. From Table 2.3.2 we see that the 
libration was computed to go to zero at 7479.1 + 35.3 = 7514.4 seconds. If the 
thruster had fired for the whole 35.3 seconds the velocity would have been 
cancelled and slightly reversed. Since the thruster must fire in the direction 
opposite to the velocity, the thrust direction must be reversed if the sign of 
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the velocity changes, This reversal of sign caused the integration to got hung 
up because the thruster started rapidly reversing direction at the point where 
the velocity went to zero. 

In order to bo;vo the problem of the thruster reversing direction, the 
program was changed to terminate the thruster firing when the magnitude of the 
libration velocity falls below a prescribed limit, The next run was done with a 
lower limit of 1 cm/sec on the libration velocity. The same initial 
out-of-plane angle of 3° was used. The thruster fired five times for 100 
seconds. The computed firing intervals for each activation were 521, 425, 345, 
234, and 136 seconds. On the sixth zero crossing the computed firing interval 
was 35.2 seconds. The thruster fired for 32.4 seconds before shutting off at 
the 1 cm/sec limit. Since the acceleration provided by the thrustor is F/m = 
2.5 newtons/550 kg = .0045 m/sec 2 , it would have required another 2.3 seconds to 
stop the motion, for a total firing time of 34.6 seconds. Successive short 
firings seem to have been either missed by the integrator or turned off by the 1 
cm/sec lower velocity limit. The run was ended at 20,000 seconds with a final 
wire length of 23.4 meters. 

The above run was repeated, lowering the velocity li.mit to .01 cm/sec. 
There were 5 firings of 100 seconds followed by a 34.56 second firing (computed 
interval 35.28 sec). The seventh scheduled firing of 3.49 seconds was skipped 
by the integrator. The eighth firing computed as 2.08 seconds was performed and 
lasted 2.03 seconds before being terminated by the 0.1 cm/sec lower velocity 
limit. The ninth scheduled firing of .31 seconds was also skipped by the 
integrator. The run ended at 20,000 seconds and the final wire length was 23.5 
meters. The out-of-plane libration had an amplitude of only 1.3 degrees at the 


end of the retrieval. 
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Starting the thruster firing at the time the system crosses the vertical 
has the disadvantage that the firing interval is not symmetrically placed around 
the time of maximum velocity. The center of the firing should be at the point 
of maximum velocity, The subroutine has therefore been changed to compute the 
approximate time of the zero crossing by adding one quarter of the libration 
period to the time of the last maximum. The firing interval can then be 
centered about the predicted time of the zero crossing. This procedure should 
work well on the last firing by leaving the systc a stopped at the zero crossing. 
If the velocity of the subsatellite is v, and the time required to stop the 
subsatellite with the thruster is dt, the satellite should travel a distance 
dtv/2 during the deceleration period. Since the thruster is fired at a time 
dt/2 before the zero crossing It should be at a distance vdt/2 from the zero 
crossing when the thruster fires and should arrive at the zero crossing with 
zero velocity. Stopping the subsatellite at a point other than the zero 
crossing does not eliminate the libration since the point at which it stops 
becomes the maximum swing angle for the subsequent oscillations, 

The next computer run was done with the same initial conditions but using 
the predicted time of the zero crossing to set the firing time. The start time 
was set equal to the predicted time of the zero crossing minus the firing 
interval of 100 seconds (the program should have used 50 seconds - half the 
firing interval) . The run was satisfactory with the exception of the error 
noted and the fact that short firings were still being skipped by the 
integrator . 

The next computer run incorporated the following changes. The start time 
was corrected to be the time of the zero crossing minus half the firing interval 
to center the firing time around the zero crossing. In order to check the 
actual time of the zero crossing, a section of code was added to compare the 
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successive values of the displacements from the vertical and print out the first 
time where a change of sign occurred (this being accurate only to within the 
stepsize used by the integrator) . The code was lati r changed to interpolate to 
find the zero crossing. The final change was to rewrite the code governing the 
start and stop times for the thruster to prevent the integrator from skipping 
over firing intervals shorter than the Integration stepsize. To do this the 
software was changed to save and make available to the thruster control 
subroutine the time of the last step taken by the Integrator. The logic was 
then split into three sections depending on whether the time of the last step 
was before, during, or after the commanded firing interval. In case 1, the last 
step is before the start time, The integrator may ask for the forces at any 
later time up to the maximum allowable stepsize, If the time is before the 
start time, the thruster force is of course zero. During the firing interval 
the thruster is on. The critical change is for times after the stop time for 
the thruster. In case 1 the program returns an "on" value for the thruster. 
This has the effect of preventing the integrator from skipping over the thruster 
firing. The integrator will try shorter stepsizes until it finds the start time 
for the thruster. In case 2, the last step is during the firing Interval. If 
the time requested is before the stop time the thruster is on, and If it Is 
after the stop time it is off. For case 3, the last step is after the stop time 
and the thruster is always off. 

The results of the computer run were as follows. Table 2.3.3 lists the 
start time, zero crossing time, end time, and firing interval (actual and 
computed) for each of the thruster firings. In the intervals before and after 
the last firing at 15,592 seconds, the out-of-plane angle was within the 
deadband of .01 radians. Since the amplitude of the out-of-plane oscillation is 
cpnstant (as measured in length units) the angular amplitude keeps increasing as 
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tho wire length decreases. Therefore, any non-zero oscillation will eventually 
exceed an angular deadband if the retrieval contiriues long enough. The stepslze 
was about 40 seconds at the time of the last firing, but the new software 
prevented the integrator from skipping over the firings. 


Table 2,3.3 

Start Time, Zero Crossing Time, End Time and Firing Interval 
for Retrieval with an Initial Out-Of-Plane Angle of 3 Degrees 
and Out-Of-Plane Thruster Firings. Times are in seconds. 


Tdeo 

Tzero 

Tend 

AT 

AT (COMPUTED) 

644 

652 

744 

100 

521 

1950 

1987 

2050 

100 

424 

3359 

3462 

3459 

100 

343 

4811 

4881 

4911 

100 

231 

6243 

6408 

6343 

100 

132 

7803 

8465 

7832 

29 

29 


9889 





11325 





12756 





14186 




15592 

16064 

15593 

1 

1 


17494 





18924 





In the runs described up to this point it has been noticed that the 
predicted times of the zero crossings do not seem to be very accurate. In the 
run just described the zero crossings listed in Table 2.3.3 are obviously not 
centered in the firing interval. The prediction is computed from the time of 
the largest amplitude which is accurate only to within the stepsize used by the 
Integrator. For the 100 sec firing, the crossings occurred 8,37,103,70 and 165 
seconds after start time, For the sixth firing of 29.4 seconds duration, the 
velocity was reduced to a small value and the next crossing was 663 seconds 
after the start time. In other words, the point where the thruster stopped 
firing became the maximum of the next libration. The seventh firing at 15592 
seconds also reduced the velocity to a low value and the next crossing was 472 
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seconds after the start time. The extrapolated time of the zero crossing 
without a thruster firing was 15616 seconds which’ is 24 seconds after the start 
time of the seventh firing. 

The thruster firings obviously alter the dynamics of the llbrations so that 
it is difficult to draw any conclusions about the natural periods from thruster 
runs. A run has been done with retrieval but no thruster to look at the periods 
of the librations. The system was given an initial out-of-plane angle of 3 
degrees and no in-plane angle. For greater resolution the simulation was run 
with output every 25 seconds rather than 100 seconds as in previous runs. The 
approximate time of the maxima were tabulated for both the out-of-plane 
displacement and the out-of-plane angle. The amplitude of the out-of-plane 
displacement was nearly constant during the retrieval , with some decrease toward 
the end as the angle became large enough to provide coupling with the radial 
variable. The time for each quarter cycle was reasonable constant at around 710 
seconds. The theoretical value for small oscillation is 694 seconds. Eor the 
out-of-plane angle the quarter cycle from the zero crossing to maximum angle was 
systematically larger (about 795 seconds) than the quai ter cycle from the 
maximum angle to the zero crossing (about 650 seconds) . The numbers were 
somewhat crudely calculated but the effect is unmistakable. There were 
fluctuations in the times for each quarter cycle of the out-of-plane 
displacement (as measured in length units) sufficient to explain the errors 
observed previously in predicting the zero crossings. Coupling with the 
in-plane component may provide an explanation but the effect was not carefully 
studied . 

In order to see just the effect of out-of-plane to in-plane coupling a run 
was done with an initial out-of-plane angle of 20 degrees and no thrusters or 
retrieval. The out-of-plane oscillation induced an in-plane libration of about 


3 degrees, A quarter cycle of the out-of-plane was around 70S seconds and there 
was a fluctuation of about 3,0 seconds, It appears that various factors such as 
thruster firings, reeling operations, and coupling between modes may complicate 
the dynamics of the tether librations so that the zero crossings are difficult 
to predict exactly. 

As an alternative to trying to predict the time of the zero crossing, the 
thrUBter control algorithm has been revised to use the tether displacement for 
setting the start time of the thruster firings. In this mode, the peak velocity 
computed from equation (4) is multiplied by half the firing time from equation 
(5) to determine at what point to start the thruster. When the tether reaches 
that position the firing start and stop times are set. A 20,000 second run has 
been done with output every 100 seconds, The tether is given an initial 
out-of-plane displacement of 3 degrees. The tension control algorithm is used 
to control the in-plane angle and the thruster to control the out-of-plane 
during retrieval. In order to better analyze the periods of the librations the 
plotting program has been revised to calculate the times of the extrema and zero 
crossings by interpolation between the output data points. For the zero 
crossing, a linear interpolation is used. For the maxima and minima a quadratic 
is fitted to the 3 points near the extremum and the time of the extremum 
computed analytically from the equation for the quadratic, The accuracy of this 
procedure depends on the interval between the output points on the plot file. 
Table 2.3.4 lists the start time, zero crossing time, end time, and firing 
interval (actual and computed) for each thruster firing. For the 100 second 
firings, the zero crossing occurs 50, 54, 57, 59 and 72 seconds after the start 
time. For the last firing of 27 seconds duration the zero crossing is 371 
seconds after the firing. This is because the motion was essentially 
neutralized by the firing. For the first firing the zero crossing is well 
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centered In the firing interval. The zero crossings occur successively later In 
the interval as the libration amplitude and velocity decrease with each thruster 
firing. 


Table 2.3.4 

Start Time, Zero Crossing Time, End Time, Firing Interval 
(Actual and Computed) , and Time of the Last Extremum for 
the Out-Of-Plane Displacement. Times are in Seconds. 


Tbeo 

Tchoss 

Tend 

AT 

AT 

(Computed) 

Thax 

604 

654 

704 

100 

521 

0 

1945 

1999 

2045 

100 

423 

1311 

3408 

3465 

3508 

100 

342 

2733 

4818 

4877 

4919 

100 

230 

4168 

6259 

6331 

6359 

100 

113 

5583 

7727 

6098 

9,331 

10959 

12389 

13619 

15249 

16679 

18109 

19539 

7754 

27 

27 

7024 


The zero crossing times in Table 2.3.4 are from the printout of the program 
and have bean determined by interpolation between steps taken by the integrator. 
These steps are shorter than the output interval of 100 seconds. In Table 2.3.5 
the results from the plotting program are tabulated. This program uses the 
output at 100 second intervals for interpolation. Table 2.3.5a is the 
out-of-plane displacement and Table 2.3.5b is the out-of-plane angle. The zero 
crossing times differ somewhat between Table 2.3.4 and Table 2.3.5 because of 
the different time intervals used for interpolation. Table 2.3.4 should be more 
accurate for this quantity. The largest difference is for the fifth zero 
crossing where the dynamics is strongly affected by the thruster firing. The 
times of the extrema given in Table 2.3.4 and Table 2.3.5a differ by a few 
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seconds. The values in Table 2.3.4 are not interpolated and are probably less 
accurate. The times for a quarter period given ’ in the last column of Tables 
2.3.5a and 2.3.5b show substantial fluctuations. In Table 2.3.5a the 
fluctuations aro only during the period of thruster firings. After the last 
firing the length of a quarter cycle is consistently close to 715 seconds. For 
the out-of-plane angle, the quarter cycles are asymmetric after the end of the 
thruster firings. It takes about 786 seconds to go from the zero crossing to 
maximum angle, but only about 643 seconds to go from the maximum angle to the 
zero crossing. Apparently the displacement is simpler to model than the angle 
during retrieval. 



Table 2.3.5a 


Tines (sec) of the Extrema and Zero Crossings 
for the Out-Of-Plane Displacement. The Last 
Column Lists the Time for Each Quarter Cycle, 


Tmh 


Tcrobb Thax Tdieb 


0 


2728 


5889 


8814 


11674 


14534 


17394 


657 

1999 

3469 

4881 

6349 

8098 

9531 

10959 

12390 

13819 

15249 

16679 

18109 


1315 

4166 

7029 

10245 

13105 

15964 

18824 


657 

658 
684 
729 
742 
697 

715 
707 
761 
680 

1069 

716 

717 
714 

714 

715 
715 
715 
715 
715 
715 
715 
715 
715 
715 
715 
715 
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Tnblo 2. 3.5b 


Times (sac) of the Extrema and Zero Crossings 
for tho Out-Of -Plano Angle. The Last Column 
Lists tho Time for Each Quarter Cycle. 

Thin Taiooa Thax Tcirr 

0 

657 657 

1352 69S 

1999 647 

2815 616 

3468 653 

4227 759 

4001 653 

5665 785 

6348 683 

7098 749 

8098 1000 

8887 788 

9530 644 

10315 785 

10958 643 

11746 788 

12389 643 

13176 786 

13019 643 

14605 786 

15248 643 

16035 787 

16679 643 

17465 787 

18109 643 

18895 786 

643 


19538 
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Figure 2.3,1 shows plots of the results for the run above, Part o) shows 
the wire length vs. time starting from 20 km and going to 23 meters, The 
subsatollito is initially stationary and slight deviations from u perfect 
exponential retrieval can be seen in the plot resulting from the action of the 
tension control law. Part b) of the figure shows the tension as a function of 
time starting from 5,6 x 10* dynes and going down to 4.5 x 10 3 dynes. The 
action of the control law is ovidont in the plot. Plot c) shows the in-plane 
displacement vs, time, Part d) shows the iri-plano angle vs. time. Under the 
action of the tension control law tho In-plano angle approaches the equilibrium 
angle of 167. 82 degrees with the amplitude of tho angular oscillations 
decreasing by about a factor of two on each cycle. The ongla of 167,02 degrees 
is 12.18 degrees from tho vortical at 100 degrees, Part o) shows tho 
out-of-plane displacement vs, time. The amplitude of the oscillation is roducad 
by the chruster firings at each zero crossing, After tho sixth firing tho 
remaining oscillations are not visible on this scale. Part f) shows the 
out-of-plane angle vs. time, Although the amplitude of tho out-of-plane 
displacement is always decreasing as seen in part e) , the out-cf-plane angle can 
increase with time as tho length of the wire decreases. The angle is still 
growing after the first four thrUBtor firings but is nearly eliminated by tho 
last two firings. Tho small remaining oscillation produces a noticeable angular 
oscillation by tho end of the run. The amplitude was within the deadband used 
of .01 radians (.57 degrees) so there were no further thruster firings, Part g) 
shows the out-of -plane vs. in-plano displacement, The lino starts at the lower 
left corner with an out-of-plane displacement but no in-plano displacement. Tho 
line ends in the bottom center of the plot with virtually no out-of-plane 
displacement and a continually decreasing in-plane displacement. Part h) shows 
the in-plano angle vs. out-of-plane angle. The line starts at the upper left 
corner with an in-plano anglo of 180 degrees (local vertical) and an out-of- 
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Eigure Captions 

Eiguro 2.3.1. Retrieval of a subsatollite deployed downward on a 20 km 
tether using tension control and out- of -plane thrusters. 

a) Wire length (cm) vs. time (sec) 

b) Tension (dynes) vs. time (sec) 

c) In-plane displacement (cm) vs. time (sec) 

d) In-plane anglo (deg) vs. time (sec) 

0) Out-of-plano displacement (cm) vs. time (sec) 

f) Out-of-plane angle (deg) vs. time (sec) 

g) In-plana displacement (cm) va . out-of-plane displacement (cm), 
0-20000 sec. 

h) In-plane angle (deg) vu. out-of-plano angle (deg), 0-20000 sec. 

1) In-plane anglo (deg) vs, out-of-plane anglo (dog), 7800-20000 sec. 
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Figure 2.3.1c 
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piano angle of "3 degrees , The in-plane angle converges to 167.82 degrees and 
the out-of-plane angle approaches zero because of the thruster firings. Part i) 
is a blow-up of the last part of Figure 2,3.1h after the last thruster firing. 
This plot covers the time period 7800 seconds to 20,000 seconds. The line 
begins in the lower center of the plot with an oscillation that is primarily in 
the in-plane direction. The sixth thruster firing has largely eliminated the 
out-of-plane oscillation. The tension control law continues to damp the 
in-plane angular oscillations, but the out-of-plane angle increases because the 
out-of-plane oscillation has a fixed amplitude (in centimeters) and the length 
of the wire is continually decreasing. At the end of the simulation the 
oscillation is primarily in the out-of-plane direction. 


2,4 Slack Tether Studies 

During the reporting period SAO has modified SLACK3 to include avoidance 
maneuvers using Shuttle thrusters and run several simulations. A number of 
approximate criteria were also developed to estimate the total possible downfall 
of tether onto the Shuttle under various circumstances. 








i. --iffrlhtrrr ,, 
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2,4.1 5LACK3: Shuttle Thruster Avoidance Maneuvers - 

The SLACKS input and setup routine was modified to accept a Shuttle 
thruster maneuver with parameters: direction {specified as pitch and roll 
angles) ; acceleration (cm/Bec*) ; initial delay (sec) and cutoff (sec) . Those 
parameters are passed in commons to BOOM, the subroutine generating the 
deployment boom tip’s position! the Shuttle center of mass position is 
computed, to which are added any rotation maneuver and boom vibration. 
Experience led us to add an option to produce plot output relative to the 
instantaneous Shuttle center of mass, rather than the coordinate origin at the 
original Shuttle c.m., passing the c.m, position in commons to the output 
routine. This is, after all, what we wish to see and otherwise the overall 
system motion dominates the tether/Shuttle motion, confusing the plot. 

We derive a Shuttle acceleration adequate for our purposes as follows*. The 
Shuttle has four forward and three reverse thrusters , with thrusts in the range 
800-900 pounds thrust; use an average of 850. Now, one "pound thrust" in 
everyday units becomes (1 pound) x (1/2 kilogram/pound) x (g = 10m/sec 2 ) = 5 
Newton. Thus one thruster produces about 850 x 5 « 4x10 3 N thrust. Taking the 
mass of the Shuttle as 10* kg, the acceleration produced per thruster is then 
4xl0 3 /10 s = 4xl0' 2 m/sec 2 = 4 cm/sec 2 . 

To examine the utility of various thruster maneuvers in avoiding the recoil 
of a broken tether we took one of the cases presented in the last quarterly 
report and added several different maneuvers. The case used was that of a 20 km 
tether cut at 1 km (see Elgure 2.3.1 of the previous quarterly). This was in 
some sense the most "dangerous" of the cases considered: high recoil velocity 
(220 cm/sec) due to full deployment, with a long remnant, (Longer remnants, 
e.g. 20 km cut at 10 km, might prove more dangerous or less amenable to 
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avoidance, duo to tho way in which tho tothor does not simply fly by the Shuttle 
but "bounces*' up and down by a fraction of n kilometer, but ns datailod in the 
lnat roport these are difficult to simulate duo to tho tether coming into full 
tension > ) 

The avoidance maneuvers all burned throe thrusters, i.o, had an 
acceleration of 12 cm/sec*. The thrust was initiated at 5 seconds after the 
break; this allows 1 (2.2 ra/sac) x (5 soc) = 11 motors of tether to downfall onto 
tho Shuttle before initiating avoidance, but same delay was rogardod as 
inevitable to avoid false alarms, Tho two other parameters each had two 
possiblo values, loading to four cases: Tho direction was either straight 
forward or 90° to one side; tho acceleration was either cut off after 20 
seconds (i.e. at t = 25) or Bimply allowed to run. The 20 sec cutoff was chosen 
to give a 20 meter displacement (about 2/3 Qhuttlo longth) at cutoff: distance 
= ijat* = 0.03t 2 meter. 

The results are shown in Figures 2.4.1 to 2.4.4. 

The two runs with no thruster cutoff terminate on excess computation time 
after having completed only a part of the requested simulation (Figures 2.4.1 
and 2.4.2) . This is probably due to the tether's trying to come into full 
tension, which as discussed in Section 2. 3. 1.3 of the last report, is a 
difficult situation for SLACK3 to cope with: it attempts many closely spaced 
bounces, and may simply fall, Effectively we are placing an artificial 
gravitational field on the tether by constantly accelerating the Shuttle, 

Overall, the runs with out-of-plane thrust are more successful in avoiding 
the tether. Thrust in this direction eliminates the complications of drag and 
Coriolis forces. Figure 2.4.3 demonstrates that a brief impulse, 20 seconds of 
thrust, is adequate for initial avoidance of the vast majority of the tether. 


0 

B 
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Elgure 2.4.1. Thruster avoidance maneuver. Original 20 km tether cut at 1 km, 
resulting in recoil velocity of 2.2 m/sec. Three thrusters, producing 
acceleration of 0.12 m/sec, initiated at 5 seconds after the break. Directed 
perpendicular to orbital plane, with no thruster cutoff. Output at 25 second 
intervals, total run of 200 seconds. Model has 25 tether mass points. 
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Figure 2.4.2. As for Fig. 2.4.1, except: 
cutoff, total run of 125 seconds. 
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1 Page 42 

I 

t (It 000 ms unlikely that we can totally avoid downfall onto the Shuttle; what is 

I desirable, and possibly achievable, is to minimize this downfall and avoid 

| t wrap-around.) However, after recoiling by the Shuttle and bouncing at full 

f 

j 

| extension in the downward direction, the tether tip "whips" around and on 

1 

! rebounding the tether comes around the other side of the Shuttle, wrapping 

* 

around it. A more complicated maneuver appears necessary to avoid this 
wrap-around (the fact that it happens in the out-of-plane direction as well as 
in-plane Indicates that we cannot simply modify the initial burn parameters) . 

Two possibilities are: (1) when the tether is fully extended below, apply 

another brief burn (in a combination of in-plane and out-of-plane directions) to 
match the Shuttle velocity with the tether tip velocity, with a bit extra for 
avoidance, and (2) rotate the Shuttle to match the average wrap-around angular 
velocity. 

2.4.2 High Resolution Loss-Of-Tension Model - 
There has been no effort on this task in the reporting period. 

2.4.3 Analytical Studies Of The Slack Tether Problem - 

I 

There has been no effort for this contract on this task during the 
reporting period. 

SAO has acquired a separate contract with Martin Marietta Corporation (MMC) 
for theoretical and innovative modeling studies of the detensioned tether. The 

majority of the analytical studies will now be performed under the MMC contract fj 

(while lumped mass studies such as SLACK3 will continue under the current NASA 
contract) . Specific applications to electrodynamic tether safety studies may be 
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time . 

2.4.4 Feasibility Arguments For Tether Impact On Shuttle - 

In the last quarterly report, Figure 2.4.2, we saw that not all severed 
tethers will recoil fully to the Shuttle, even ignoring the Coriolis and drag 
forces which displace them from the direct path. In the particular case 
illustrated there, a 20 km tether, cut at 10 km merely "bounced" toward the 
Shuttle and away by a half kilometer or so. 

In this section, we present some crude, "back-of-the-envelope" arguments 
addressed to questions such as: Will the tether go slack at all? How much of 

the tether might conceivably impact on the Shuttle? 

Throughout, we assume that the tether is initially deployed straight up (or 

down) and is in equilibrium. The following notation is used: 

t - original tether length 

r = remnant tether length (after cut) 

D = downfall tether length (possible impact on Shuttle) 

m = subsatellite mass 

/i = tether mass per unit length 

EA = tether elasticity 

c = tether speed of sound (= \/EA/n ) 

Cl = orbital angular velocity 

In general we make the crudest possible assumptions; e.g. we use only the 
subsatelllte mass, and not the tether mass, when computing the gravity gradient 
force before the break. We also assume throughout that the amount of tether 
stretch is small compared to the natural length. In numerically evaluating the 
formulae, we use Cl ~ 1.2 x 10" 3 sec" 1 , ra - 500 kg, fi = 1.1 x 10" 3 kg/m, EA = 0.6 


x 10 s kg m/sec z (c = 2.34 x 10 3 m/sec) . 




T 
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For our first criterion* we ask the question "Is there enough stored 
elastic energy in the tether remnant at the moment of the cut for the entire 
tether to move far enough against the gravity gradient field to become 
detensioned?" , The significance of this question is based on the observation 
that the wave of slackening proceeds along the tether from the point at which it 
was cut, as seen in Section 2.2.4 of Quarterly Report #2 with the high 
resolution loss of tension model. That simulation did not include gravity 
gradient forces , but in their presence we should expect that the sharp tension 
profile would become smoother as the wave propagates toward the tether j as the 
gravity gradient potential energy dominates the elastic energy, the detensioning 
wave will at some point merely be reducing the tension to some positive value 
(rather than zero) by the time it reaches the Shuttle, 


To quantify this argument we need two energies: The stored elastic energy 

in the remnant, Eel- The energy, E B l, needed to move against the gravity 

gradient force from the initial configuration to one in which each element of 

the tether is just barely detensioned, For Eel, we assume a massless tether; 

this simplifies the tension calculation and also means that each tether element 

is stretched by the same amount. The tether spring constant of the segment 

which will become the remnant is then k = EA/r and the gravity gradient forcB is 

. . „ „ 3fi 2 m£r 

E — 3n 2 fm. The total stretch of the remnant portion is then Ar = F/k ~ — — — . 

Considering this segment as a simple spring, the stored elastic energy is then 
9 n 4 m 2 


Eel - hk(Ar ) 2 = - 


rf z 


To compute Esl, we first compute the energy 


2 EA 

required to raise each mass element of the tether from its position in the 
initial configuration to that in the "barely slack" configuration. Let s be the 
natural (detensioned) length along the tether from the Shuttle attachment to a 
given element. The tether has been uniformly stretched by an amount e = A r/r = 
3ft z m£/EA. Consider an element of length ds and mass dm = (ids at coordinate s. 
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It was originally displaced by an amount es, and now must move by this amount 
against the gravity gradient force, which is F (s) - 3n 2 sdm. The energy required 
is thus dE - (es)F(s) = 3«0*/is 2 ds, where we have used dm = /ids. The total 
energy required to move from the initial configuration to the zero-tension 
configuration is then obtained by integrating dE from s=0 to s=r and 
substituting for e: Ebl = 3 r3 Note the r 3 dependence of Eat, (in a 

uniform gravity field it would be r 2 ) and the simple r l dependence of E Et ,j this 
implies that for a given original length l, the energy required to go completely 
slack will be negligible compared to the stored elastic energy for short 
remnants, but as the remnant length is Increased the elastic energy will at some 
point be unable to overcome the gravity gradient forces. This cross-over point 
is found by equating the two energies : 

r - {W l (1) 

= >/75km l 


Eor shorter remnant lengths, the tether has enough elastic energy to go 
completely slack; for longer remnants, there is not enough elastic energy to go 
completely slack. (A short remnant, however, is not guaranteed to go slack.) 
As discussed above, if a remnant does not go completely slack, it seems likely 
that the portion near the Shuttle will remain taut. Note that any remnant from 
a tether with original length less than 75 km will have the energy to go slack, 
since the critical remnant length will be longer than the original length. 

Our next two criteria respond to the question "Assuming the tether becomes 
slack, how much of the tether could conceivably impact on the Shuttle?" In both 
cases we Ignore drag and are only concerned with whether or not such Impact is, 
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In some sense, energetically allowable. 

For criterion number two, we imagine that at the moment the tether goes 
completely slack it is cut up into infinitesimal mass elements, These elements 
will all be moving toward the Shuttle with a uniform velocity given by ^rv J = 
Eel i leading to v = 3fl 2 mf/^c where c s \/EA /(jl , In a Shuttle-centered 
coordinate system with the x-axis pointing radially outward and the y-axis 
forward along orbit, the governing equations for a given mass are: 

x as 3tl*x + 2fly 
y = - 2nx 

which can be solved (see Section 4.2.1 of the Einal Report for NAS8-35036, 
"Investigation of Electrodynamic Stabilization and Control of Long Orbiting 
Tethers, 1 ' 1 PI G, Colombo) to give 

x (r) 53 s (4 - 3cos(r)) - Vsln(r) 
y(r) = 2V (1 - cos(t)) + 6s (sin(r) — r) 

for b. mass element with initial conditions x(0) = s, y(0) =0, x(0) = -v = -OV, 
y (0) = 0 j r is a dimensionless time, r = fit. 

We are interested in which elements cross the Shuttle's orbit, i.e. have 
x(r) < 0 for some r. The maximal s for given recoil velocity V which allows 
this is the downfall D. x(r) = 0 can be written 

V 

4 - 3cos(r) = — sin(r) 
s 


and by finding the value of V/s for which the curve on the left is just tangent 
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w 




to tho curve on tho right we achieve the downfall D. Thia is 


V 

D 


vr 


which gives rise to 


D = 


30m 


V e 


t = 0.026 £ 


( 2 ) 


Note that this estimate of the downfall is independent of where the tether is 
cut. 


For the third criterion, we assume that the elastic energy E&, is much 
greater than needed to achieve full slackness. Once slack, we assume that the 
tether moves toward the Shuttle uniformly, and that any tether Impacting on the 
Shuttle is retained. We neglect drag and Coriolis forces, so that the tether 
continues to fall directly on the Shuttle. How far can the entire tether move 
toward the Shuttle against the gravity gradient force before exhausting the 
stored elastic energy which has been converted to kinetic energy? The potential 
energy of a tether of length £ in the gravity gradient field is Eco(L) = 
l 3 , We then wish to find the downfall D from a tether cut at length r 

such that Eel - E^ (r-D) - Eco (r) . If we assume that the downfall is only a 
small part of the remnant length, then we have 


D = 3 


n 2 m 2 £ 2 


/iEA 


(3) 


= 1.6 x 10 


-3 


£ 2 


Now this criterion, in contrast to the free orbit criterion (2), depends on both 
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the original length t and tho romnant longth r. Wo can put two limits on this: 
first, lot the tether be cut at its full length, r = £, which gives directly: 

D = 0.0016 t 

and second, after noting that the downfall increases as tho romnant longth 
decreasos , suppose that the entire remnant becomes downfall, leading to D a - 
0.0016 £2 or 


D - 0.040 t 

This latter estimate agrees well with that of (2) numerically, and indeed is 
only different by a factor of s/ 3/7 . 

How do these compare with the behavior observed in the last report? There, 
a 20 km tether cut at 0.5 and 1.0 km recoiled fully, When cut at 10 km, it only 
bounced up and down by about ij km or so. For a 20 km tether, tho various 
estimates of downfall are 0.6 and 0.8 km for complete downfall [free orbit (2) 
and D = r in (3)] and 0.07 for a long remnant [r ~ 10 km in (3)] . Our estimates 
appear somewhat conservative, though not excessively so. They could probably be 
somewhat refined by different or more sophisticated arguments. Perhaps, for 
instance, one could include "flyby" as well as "downfall": much of the tether 

actually bypasses the Shuttle, and as it falls down the gravity gradient 
potential on the other side it effectively pulls the remaining remnant behind it 
rather than simply vanishing from the equations as in the derivation of (3) . 
Also, it would be Interesting to conduct a series of simulations for remnants 
between 1 and 10 km to determine where the difference in behavior occurs. 
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2.4.5 Concluding Remarks - 

The slack tether simulation program SLACK3 has been extended to include 
Shuttle thruster maneuvers. Simulations have boon made Incorporating thrusters 
for avoidance maneuvers. Simple brief bursts perpendicular to the orbit 'plane 
provide good initial avoidance, but prevention of long-term wrap-around will 
require more complex maneuvers. Design and simulation of those long-term 
avoidance maneuvers is an appropriate task for the next reporting period, 

Several criteria for the impact of tether remnants on the Shuttle, based on 
energy arguments and similar general considerations, have been derived. The 
numerical results appear somewhat conservative when compared with the 
simulations reported in the previous quarterly report. The arguments leading to 
the criteria should be refined, and further simulations should be run to gain 
more precise knowledge of the transition from the complete downfall/flyby of 
short remnants to the bouncing behavior of long remnants. 

Detailed theoretical (non-ball-and-spring) studies of the slack tether are 
now being carried out under separate contract. Appropriate results of these 
studies may be applied to specific cases of Interest to the current contract. 
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2.5 Numerical Calculation Of The Electric Eiold Around An Electrodynamic Tether 
After Breakage 

2.5.1 General - 

In Quarterly Report we deecribed the development of a numerical program 
which solves Laplace's equation by the Liebmann net method. It was explained 
there that a geometrically varying grid was used to span a large region of space 
while allowing for a fine mesh spacing in the neighborhood of the end of the 
tether. The potential at the boundary of the grid was initialized by an 
analytical calculation of the field of a metallic prolate spheroid. This 

potential is expressed in terras of the Legendre function of the second kind of 
order one, Qi(£) whore, as explained in Quarterly Report 

Qi-€Qo(€)-i (l) 

and 

Qo (O = m (2) 

The prolate spheroidal coordinate £ is defined as 

< = ri l^' < 3) 

The distances from the two focal points of the spheroid to any point in space, 
ri and rj, are given in terms of the cartesian coordinates (x,y,z) by 


n = >/x z + y z + (z - c) 2 


(4a) 
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ra « >/x* + y* + (z + c) 1 (4b) 

In this analysis as a result of axial symmetry wo have replaced x and y by r , 
where: 

r = v'x* + V* (5) 


Thus 


n => s/r* + 

(*-<=)* 

(6a) 

r 2 = s/r* + 

(Z + c) 2 

(6b) 


The numerical instability we reported in Quarterly Report is a direct result 
of the square root operation. When r is small with respect to ( 2 +c) or (z-c) 
then large round-off errors occur, and when z approaches c for large r equation 
(6a) causes large round-off errors. 

When one term is much smaller than the other one can use the approximation 

y/I + A at 1 + | (7) 

which is easily tested by squaring both sides of (7). If we substitute (6a) and 
(6b) into (2) we find that the argument of the log function becomes 

ri + r a + 2c 
ri + r 3 - 2c 

It Is essential that care be taken to distinguish between the two cases z > c 
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f 

I and z < c when using the approximation (7) with (8) . Since rj, and rj are 

i ' 

positive numbers, in the limit as r -» o, \/ (z~c) * -» |z-c |. When jz-c j is 

i 

f ( smaller than r wo find thnt (7) results in 

(9) 

When A in approximation (7) is not sufficiently small the error Introduced by 
the approximation is unacceptably largo. We have reisolvod this difficulty by 
replacing \fx +~x with its Taylor series expansion. 

If; 


ri a I z - c I + -t -~ — r 
' ' 2j z - c ] 


f = \/l + x 


(10) 


then : 


f 


1 + 


x 

2 


x z , 3 x 3 

8 + 8 31 


15 x< 

16 4T + 


( 11 ) 


This procedure accurately calculates \/l + A even when A -+ In the following 

computation up to 30 terms in the series are used. This technique has 
completely removed the numerical instability from our computer model . 
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2.S.2 Computer Model Of The Eloctrodynamic Tether - 

2 .5. 2.1 Grid Calculations - 

The mesh box width A(n) are increased in a geometric progression such that 

A(n) « a A(n - 1) (12) 

<* 

Thus n mesh boxes span a distance t whero 

t = A(l) [l + a + a 2 + a 3 + , , , a"* 1 ] ( 13 ) 

£ " 4<1 > ?=T (14) 

To specify t and n and calculate a would retire the solution of a high order 
polynomial, We devised a simple numerical algorithm to calculate a given l and 
n} such algorithm is implemented in this software code, 

2. 5. 2. 2 Contour Plots - 

The grid chosen to model the tether covers a region ± 2 kilometers from the 
end of the tether in the z direction with a radius of 2000,001 meters. It is 

* 

defined as an array of 91 x 161 points for a total of 14,400 mesh boxes. An 
equipotential contour plot of the analytically calculated field of a conducting 
prolate spheroid which models approximately the tether embedded in an axial 
field of 0.169 volt/m is shown in Figure 2,5.1. The region of the grid occupied 
by the tether was then set to zero and an index array which inhibits the 
calculation of potentials on the tether was initialized. A contour plot of the 
potentials calculated numerically by the Liebmann net method for the cylindrical 

| 
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tothor (1 mm radius and 20 km long) are plotted Cor the outermost 2 km of tho 
tether in Figure 2,5.2. These two calculations are compared by superposing the 
two contour plots in Figure 2.5.3. Because of the nature of the high ratio 
geometric grid used for this calculation the shape of the field at the end of 
tho tether is obscured in these plots. To resolve this difficulty the plotB 
algorithm was modified to "zoom" in on the end of the tether. The first stage 
grid consisted of 90 X 160 mesh boxes, the second stage 50 x 80, the third stage 
30 x 40 and the fourth 20 x 20. Thus the axial lengths of the grid in these 
four stages are 4 km, 1.9 m, 40 mm and finally 5 mm. Stages 2 through 4 are 
shown in Figure 2.5.4 through Figure 2.5,6. 

2.5.3 Numerical Results - 

This calculation was performed for the 28° orbit in which the 20 km tether, 
2 mm in diameter, is exposed to a magnetically induced field of 0.189 v/meter 
(the worst case) . The result of the calculation (vacuum case) are best shown in 
Figure 2,5.6. The voltage drop across the Insulator, 2.5 mm from the tip of the 
tether (point B) , is about 100 volts while the radial and the axial component of 
the electric field are respectively given by: 

E r = 330 Kvolt/m (15) 


E r = 0 

At the corner of the tether tip (point A) , where we have assumed a squared edge 
and a flat top of the wire, the electric field is: 


E r = 1139 Kvolt/m 


(16) 










SUPERPOSED POTENTIAL PLO 
A V: 200.0 ZSPAN: ML . 0 


Figure 2.5.3. Comparison between analytically ai 
potentials around a 20 km long tether, 
superimposing Figure 2.5.1 and Figure 2.5.2. 
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E, m 1350 Kvolt/m 

These are bIbo the maximum values for the electric field in the case under 
investigation. At the very tip of the tether (point C) , whore the tether is 
assumed to be flat, the electric field is: 

E r = 0 (17) 

Ex = 700 Kvolt/m 

This computation takes into account the effect produced by the Insulator around 
the conducting wire. The insulator characteristics for the TSS first mission, 
provided by MMA, are as follows t 

insulator thickness ~ 0.36 mm 

dielectric constant = 2.1 (18) 

dielectric strength = 1200 volt/mil = 47 Kvolt/mm = 47000 Kvolt/m 

It is therefore interesting to notice that (in vacuo) the dielectric strength of 
the insulator is never exceeded even after the tether breakage. It is also 
interesting to see how the value of the electric field computed by this high 
fidelity model compares to the simplified single-dimensional computation 
performed in Quarterly Report #1. In that Quarterly Report we prellminarly 
computed the electric field at point C in Figure 2.5.6 assuming that the tether 
was a prolate spheroid. The prolate spheroid however has a too small radius of 
curvature at the tip so that we corrected that value by a factor of 2 x 10 7 that 
relates the radius of curvature of the prolate spheroid to that of a tether with 
a semispherlcal tip. The result was as fallows: 




w 

r 
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E. « 4.5 x 10 l3 /20 7 = 225 K'.volt/m 


(19) 


This value i'« of the same order of magnitude of that computed with the high 
fidelity model. 


2.6 Preliminary Estimation Of The Electrodynamic Hazards Due To The Breakage Of 
A Tether Embedded In A Plasma 

2.6.1 Introductory Remarks - 

A preliminary analysis is underway of the mechanisms by which a 
Shuttle-borne long elsctrodynamic tether, broken at a certain point along its 
deployed length, may produce harmful effects. 

2.6.2 Effects Of Break - 

Assuming a break occurs near the Shuttle, two of the possible hazardous 
effects are the following: 


\ 

i 

i 

I 

* 

.fii 


l! 

1 


r 


a) Appearance of high voltage across the break, which can accelerate 
electrons toward the Shuttle, so that the Shuttle is bombarded by high-energy 
electrons , 

b) Arc drawn in vacuum between separating contacts extinguishes abruptly, 
radiating a pulse of electromagnetic energy that may be powerful enough to 
produce EMP"type effects. 

We consider, first, effect (a), then (b) . 
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2,6,3 High Voltage Across Break - 

When a circuit carrying current is interrupted, either intentionally by a 
switch or inadvertently by a break as in the present case, an arc is drawn 
between the separating contacts (Cobine, 1958) . The present problem represents 
the caso of a vacuum arc (Lafferty, I960) . As the contacts separate and the gap 
length increases, a point is reached where the available voltage is insufficient 
to sustain the arc and the arc is extinguished. If there were no Inductance in 
the external circuit, then at this gap length the post-extinction voltage across 
the gap would rise from its normal operating value (vicinity of 20V) to the 
circuit driving voltage or e.m.f. (Cobine, 1958), which is about 4 kV in the 
case of the ttlectrodynamic tether, 

If there is inductance in the circuit, then the current continues to flow 
while the contacts continue to separate further. At some (presently unknown) 
point the arc extinguishes, and the voltage across the gap simultaneously rises 
to its "extinguishing" value, which depends on the arc current-voltage 
characteristic. The extinguishing voltage can be much larger than the circuit 
driving voltage (Cobine, 1958) , 

In addition, because of stray capacitance in the vicinity of the arc, the 
interruption can give rise to high transient voltages across the contacts. This 
is because the magnetic energy stored (LI z /.2, where L is the inductance and I is 
the current) flows into the capacitance in the form of electric energy (CV 2 /2, 
where C is the capacitance and V is the resulting transient voltage) . This 
energy subsequently oscillates between the two forms. If all the energy is 
transferred in this way (no losses) then, equating the two energies, the 
resulting voltage across the capacitance (and therefore across the contacts) 


would be: 
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V » I (L/C) 1 /* 

where (L/C) 1 /* is the "surge impedance," and I is the current flowing prior to 
interruption. Thus, the voltage can be large independently of the driving 
voltage If the surge impedance is large. There is also a resonant frequency f 
determined by the relation: 


2irf » (LC)'V* 

Hence a small stray capacitance results in a largo surge voltage, and an 
oscillation at high frequency. 

Wa .oed therefore to address the problem of the evaluation of how much 
voltage can appear across the, contacts and thus give rise to accelerated 
electrons that can bombard the Shuttle. Some rough numerical estimates will be 
given later. 

2.6.4 Pulse Of Electromagnetic Energy - 

As the contact separation increases, the arc (drawn in vacuum) grows in 
length. At some point (not presently known) the arc extinguishes. The 
interruption of current results in electromagnetic radiation, the properties of 
which depend on the details of the time and space behavior of the current. 
Rough estimates of this radiation pulse are given in the next section. 


lb. 
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2,6,5 Numerical Estimates Of Bombardment - 

Assume the following parameters i 
Tether length ** 2 x 10 4 m 
Conducting wire radius; 0,5 mm (rj) 

Outer radius of wire insulation: 0.75 mm (rj) 

Inductance of 20 km tether: 0.012 h (L) 

Capacitance of 20 km wire (with respect to rj) : 2.0 x 10' 6 F 

Tether current: 1A 

Break separation velocity: 4.2 m/s 

Electron conical beam semi-angle: 30 degrees 

We may obtain estimates of the surge impedance (L/C) and therefore of gap 
voltage if we know the values of C and L. Since the problem is a dynamic one 
involving not only the arc itself but also the tether with finite propagation 
times and attenuations of waves moving along it, the appropriate effective 
values of L and C are difficult to estimate. One extreme limit is obtained by 
treating the whole tether as a lumped circuit element with L = 0.012 h and C — 
2.8 x 10' 6 F. This gives an underestimate, namely, a surge impedance of 65 £1, 
and therefore a low surge voltage of only 65 V. 

At the other extreme, an overestimate is afforded by taking the capacitance 
as that of the wire tip, namely, 5.6 x 10' 1 * F, while retaining the whole tether 
inductance, 0.012 h. This gives a surge impedance and voltage of 4,6 x 10 s 0 
and 460 kV. While admittedly this overestimate of V results from using too 
large a value for inductance, namely, the whole-tether inductance, it is useful 
as a worst-case since the appropriate value of inductance is unknown. A 
reasonable compromise value for V may be obtained as the geometric mean between 
the lower limit, 65 V, and the upper limit, 460 kV, namely, 5500 V. (Note that 
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thin value is comparablo with the normal driving voltage, 4000 V.) 

2.6.6 Bombardment Of The Shuttle - 

The time during which electrons would bo accelerated at this voltage is 
given by 1/2 of the period of the high-frequency oscillation , namely, jr(LC) l /* ( 
or BO ns, During this time period, 5-kov electrons, moving with velocity 4.2 x 
IQ 9 cm/s can travel 3,4 m, However, henceforth wo will assume a distance of the 
order of 1 meter. 

To estimate the particle energy flux onto the Shuttle, wo assume that 1 
ampere of electrons, confined to a conical beam of semi-angle 30°, strikes the 
Shuttle one meter away. The angle 30° is characteristic of beams emitted from 
the cathode spots in vacuum arcs (Reece, 1963) , Assuming that one ampere of 
5-kev electrons hits one square meter of Shuttle surface, we obtain 6.25 x 10 14 
eloctrons/cm*/s • This is orders of magnitude larger than the fluxes of 
radiation-belt electrons of comparable energy. Hence, some damage may be done 
by the pulse, but this needs yet to be quantified. 

2.6.7 Numerical Estimates Of EMP - 

The electromagnetic radiation depends on the details of the time and space 
behavior of the current. In particular, the current distribution in the arc is 
of interest, but this is an unsolved problom to date (see e.g. Lafferty, I960). 
Eor example, a conical jet of metallic vapor is emitted from a cathode spot. 
Electrons are also emitted by thermoionic field emission. Avalanches occur, 
producing ionization of the metal. Some of the ions, but not all, bombard the 
cathode releasing more electrons. The dense metallic vapor moves at high speed, 
dragging along the ions. At some point, the jot has expanded until it is so 
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rarefied that the ion mean free pnths are long and the ions are no longer 
drag^iyi along. Thus, the plasma is highly collision-dominated near the cathode, 
and becomes collisionless downstream. This transition between 
collision-dominated and collisionless regimes is partly responsible for the 
difficulties of modeling this phenomenon, 

However, we will make a rough estimate, as follows. Assume that the arc 
current is one ampere, and that it is one motor in length before it 
extinguishes. Assume that the extinction of current 1 occurs over a period of 
10-100 ns (Earrall , in Lafforty, 1980, p. 194-5). We will use a dipole 
approximation. Assume the length of current is h - lm, The radiation electric 
field at n distance of r may be approximated by 


E = b dl/dt 
47f<C s r 

If I changes by 1 A in 10 ns, dl/dt "10 s A/s. Using this, together with h 
= lm, e = 8.85 x 10‘ 13 E/m, C <= 3 x 10 s m/s, we obtain 

E a 10/r (V/m) 

Hence , at r - 1 m the induced electric field is about 10 V/m, l.e,, a weak 
field. This may be an indication that the EMP-type threat is weak. However, 
the model is too crude, and needs to be examined more carefully, 

It should be noted that the energetic electrons discussed above may 
liberate x-rays due to the surface bombardment. These x-rays may possibly be a 
hazard that warrants further examination. 
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3.0 FROBLEMS ENCOUNTERED DURING REPORTING PERIOD 




None , 


4.0 \CTIVITY PLANNED EOR THE NEXT REPORTING PERIOD 

The simulation activity on the satellite's rotational dynamics will 
continue by performing retrieval simulations with and without out-of-plane 
tether libration and thruster libratlon control. An investigation of the 
coupling between thrusters' activation and tether dynamics, in the case of 
constant tether length, will be Initiated. 


The analysis of Shuttle's maneuvers in order to avoid the recalling tether 
after its severance will be further developed. Simulation of most important 
cases, as agreed upon by NASA/MSEC experts, will be carried out. Optional 
solutions to minimize dangers to the Orbiter will be looked for. 


The study of electrodynamic hazards caused by the breakage of the TSS 
electrodynamic tether embedded in a plasma will continue by improving on the 
understanding of the most likely consequences. 









